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Abstract. We use the COBE/DIRBE (1.2, 2.2, 60, 100, 140, and 240 fim) maps and the COBE/FIRAS spectra (for 
the wavelength range 100 - 1000 fim) to constrain a model for the spatial distribution of the dust, the stars, and 
the gas in the Milky Way. By assuming exponential axisymmetric distributions for the dust and the stars and by 
performing the corresponding radiative transfer calculations we closely (given the simple geometry of the model) 
reproduce the FIR and NIR maps of the Milky Way. Similar distributions for the atomic and molecular hydrogen 
in the disk are used (with an inner cut-off radius for the atomic hydrogen) to fit the gas data. The star formation 
rate as a function of the Galactic radius is derived from the FIR emission and is well in agreement with existing 
estimates from various star formation tracers. The gas surface density is plotted against the star formation rate 
density and an "intrinsic" Galactic Schmidt law is derived with excellent agreement with the "external" Schmidt 
law found for spiral galaxies. The Milky Way is found to consume ~ 1% and ~ 10% of its gas in the outer and 
inner regions respectively (for a period of 0.1 Gyr) to make stars. The dust-induced B-V color excess observed in 
various directions and distances (up to ~ 6.5 kpc) with well-studied Cepheid stars is compared with the model 
predictions showing a good agreement. The simple assumption of exponential distributions of stars and dust in 
the Galaxy is found to be quite instructive and adequate in modeling all the available data sets from 0.45 /im 
(B-band) to 1000 /xm. 
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J> ■ 1. Introduction thus rendering their use quite limited and uncertain (see, 



however, Kewley et al. 2002 for a discussion on properly 
Our Galaxy, the Milky Way, constitutes the best labora- accounting for extinction). The far infrared (FIR) emis- 



. tory for studying the properties of the interstellar medium ^ if modeled proper l y to distinguish between the dif- 

(ISM) m spiral galaxies. Being inside the Galaxy, the ob- fuge dugt emission and the emission coming from the star 

server has the advantage of looking through different lines formation complexes, has been proven to be one of the 

of sight passing both through dense environments (e.g., mogt promismg SFR indica tors (see Misiriotis et al. TOI 

the Galactic center) and through areas almost free of dust and for & review Kylafis & Misiriotis ggrj^ . 
and gas (e.g., the Galactic poles). The way that dust and 

gas (both in molecular and neutral phases) are distributed Some efforts have been made in modeling the interstel- 

inside the Galaxy not only shapes the appearance of the lar medium (ISM) in the Milky Way. Kent et al. i|1991[> 

Galaxy in different wavelengths, but also drives star for- used the K-band data from the SPACELAB infrared tele- 

mation, a key process that determines Galactic evolution, scope and modeled the distribution of the light in the 

There are several tracers for measuring the star for- Galaxy by assuming an exponential distribution for the 

mation activity in a galaxy (e.g., ultra-violet (UV) cmis- starlight in the disk and the bulge. Using the COBE 140 

sion, far- infrared (FIR) emission, H Q , and radio emission), and 240 /xm maps along with radio surveys of HI and 

yet it is far from trivial to get an accurate estimate (see H 2 , Sodroski et al. I|1994|) modeled the FIR emission of 

Kennicutt. I1998bl for a review). Most of the star forma- the Galaxy and derived the dust temperature profiles as 

tion rate (SFR) tracers are affected by dust extinction, well as gas-to-dust mass ratios and emissivities along dif- 

ferent lines of sight. A composite of COBE/DIRBE and 

Send offprint requests to: Angelos Misiriotis, e-mail: IRAS/ISSA maps made by Schlegel et al. fEHB} provides 

angmis@physics.uoc.gr a unique tool for estimating the Galactic extinction along 
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different directions. Davies et al. (|1997ll modeled the FIR 
emission from the Galaxy by fitting exponential dust dis- 
tributions to the COBE/DIRBE maps. In this study, the 
authors find that the best fit could be achieved by assum- 
ing that the dust distribution is more extended than that 
of the stars, namely that the scalelength of the dust is 1.5 
times that of the stars and the scaleheight of the dust is 
twice that of the stars. Modeling of nearby edge-on galax- 
ies (e.g., Xilouris et al. |1999t and references therein) do 
support the finding of the dust having a larger scalelength 
with respect to the stellar one, but not the scaleheight 
relation found in this study. In our analysis we will per- 
form modeling of the near-infrared (NIR) and the FIR 
maps of the Milky Way and thus we will investigate the 
relation between the dust and stellar geometrical char- 
acteristics. Drimmel & Spergel 1)2001[) presented a very 
detailed modeling of the Milky Way dust content. Based 
on a total of 48 parameters (26 parameters for the dust 
and 22 parameters for the stars), this model quite nicely 
fitted the NIR (J-band) and FIR (240 /im) radial profiles. 
In a subsequent paper (Drimmel et al. I2003[) . this model 
was used to calculate the extinction to any point within 
the Galactic disk. This model provides a detailed descrip- 
tion of the Galactic morphology, taking into account the 
structures of the disk, the spiral arms, and the warp. In 
contrast, our approach aims to model the Galaxy using as 
few parameters as possible. We know in advance that we 
will not be able to get an excellent representation of the 
detailed morphology, but we will be able to draw results 
on the large-scale distribution of the dust, the stars, and 
the star formation in our Galaxy. 



In this paper we use a simple (with the minimal 
number of parameters required), but realistic, three- 
dimensional model, which includes the effects of absorp- 
tion and scattering of the stellar light by the dust in 
the optical and near infrared (NIR) wavelengths. We also 
take into account the emission from the diffuse dust and 
the emission from a warm dust component associated 
with the star-forming regions. In order to constrain the 
model parameters, we fit the model to all the available 
data simultaneously. The data that we use consist of the 
COBE/DIRBE (1.2, 2.2, 60, 100, 140, and 240 /im) maps 
and the COBE/FIRAS spectra (for the wavelength cover- 
age between 100 and 10000 /im). Using published maps for 
the atomic and molecular hydrogen, we model the gas dis- 
tribution. After calculating the SFR density throughout 
the Galactic disk and comparing it with the observed gas 
surface density, we investigate the validity of the Schmidt 
( 1959} law as quantified by Kennicutt l|1998af) . Finally, us- 
ing a catalog of Cepheid stars with good estimates of their 
distance and B-V color excess, we confirm the validity of 
our model in the optical region. 



2. The data 

2.1. The COBEdata 

The Cosmic Background Explorer (COBE) satellite, 
launched on November 18, 1989, has been proved to be a 
very successful space mission (Boggess et al. I1992f) . With 
the Diffuse Infrared Background Experiment (DIRBE) in- 
strument, COBE surveyed the entire sky in 10 photo- 
metric bands covering the wavelength region from 1 to 
240 /im at an angular resolution of 0.7 square degrees 
(Sodroski et al. 1994). Using another instrument, the Far 
Infrared Absolute Spectrometer (FIRAS), COBE mea- 
sured the spectrum of the dust emission of our Galaxy in 
the wavelength range from 120 (im to 10000 /im (Fixsen et 
al. H996[) . These two instruments provide us with a unique 
set of data of the dust emission for a wide range of wave- 
lengths, ideal for studying the properties and the distri- 
bution of dust in the Galaxy. 

2.1.1. NIR and FIR maps 

DIRBE was designed primarily to conduct a systematic 
search for an isotropic cosmic infrared background in 10 
photometric bands from 1.25 to 240 /im (Boggess et al. 

Hauser et al. lT9T)8) . Boggess et al. l(T9"§!S|) report an 
rms sensitivity per field of view (0.7 square degrees) of 1, 
0.9, 0.4, 0.1, 11, and 4 (10 9 Wm^si-- 1 ) at 1.2, 2.2, 60, 
100, 140, and 240 /im, respectively. 

For the purposes of this study we make use of the 1.2, 

2.2, 60, 100, 140, and 240 /im Zodi-Subtracted Mission 
Average maps obtained from the COBE database. This 
set of data consists of weekly averaged intensity maps with 
the zodiac light subtracted from the maps (see Kelsall et 
al. 119981 for a description of the model for the interplan- 
etary dust). The data were extracted from the original 
data-cubes and analyzed using the special COBE analysis 
software UIDL 1 written in IDL. 

2.1.2. FIR and submillimeter spectra 

FIRAS is a polarizing Michelson interferometer (Mather 
1982) with two separate spectral channels (the high- 
frequency channel, extending from 120 to 500 /im and 
the low-frequency channel extending from 500 to 10000 
/tm). The in-orbit absolute calibration of FIRAS was 
accomplished by using an external black-body calibra- 
tor. The spectra, already corrected for Cosmic Microwave 
Background and zodiac light contribution, were extracted 
from the original data-cubes and averaged for eight dif- 
ferent directions on the sky over a range of 20 degrees in 
Galactic latitude (from -10 to 10 degrees) and from to 
10, 5 to 15, 15 to 25, 25 to 35, 55 to 65, 85 to 95, 115 to 
125, and 135 to 145 degrees in Galactic longitude. 



http:/ /lambda. gsfc.nasa.gov/product/cobe/cgis.cfm 
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2.2. The HI and H 2 maps 



3.1.1. The warm dust component 



The Leiden/ Argentine/Bonn survey of the Galactic HI 
(Kalberla et al. I2005fl is a merging of the Argentino de 
Radioastronomia and the Leiden/Dwingeloo surveys. The 
angular resolution of the combined survey is 0.6 square 
degrees (comparable to that of the COBE maps) and is 
the most sensitive, to date, HI map of the Milky Way. 

A composite CO survey of the Milky Way constructed 
from 37 individual surveys has been carried out by Dame 
et al. (2001) 2 . The total area covered by this survey is 
9,353 square degrees, which accounts for nearly one half 
of the area within 30 degrees of the Galactic plane. Using 
a CO-to-H 2 mass conversion factor of A" = 1.8 ± 0.3 x 
10 20 cm" 2 K" 1 km" 1 s (Dame et al. EODT) . a H 2 map is 
constructed. 



3. The model 

The model that we use quite accurately simulates the 
physical processes (see Kylafis & Bahcall 119871 and 
Popescu et al. I2000f) that take place inside a galaxy for 
the wavelength range between optical and submillimeter 
(submm). In the optical/NIR wavelength range, absorp- 
tion and scattering of the stellar light by the dust grains 
occurs, while emission from the heated dust grains is tak- 
ing place in the FIR/submm wavelengths. The different 
nature of the processes allows us to treat the two wave- 
length ranges separately as we describe below. Since the 
dust reveals itself by its FIR/submm emission, it is reason- 
able to begin with the modeling in these wavelengths to 
constrain the dust parameters. By fixing the dust param- 
eters in the radiative transfer model, we can then model 
the optical/NIR emission of the Galaxy. 

For computational reasons, all the radiative transfer 
calculations are made inside a cylinder of a radius of three 
times the scalelength of the dust and a half-height of six 
times the scaleheight of the stars. 



3.1. The FIR/submm wavelength range 

The emission model that is used for the FIR/submm wave- 
lengths is along the lines of the model described in Popescu 
et al. ( 2000 see also Misiriotis et al. 2001 ). Following these 
studies, we assume that the Galactic dust can be described 
by two components. The first component is the warm dust 
associated with star formation complexes, which is heated 
locally by the young stars. The second component is the 
cold diffuse dust, which is heated by the diffuse radiation 
field in the Galaxy. For the dust grain cmissivity in these 
wavelengths, we use the values reported in Weingartner & 
Draine pUUTj) . 



For the density of the warm dust component we assume 
that its three-dimensional distribution is given by the fol- 
lowing equation 



i R \ z \ 

p w {R, z) = p w {0, 0) exp ( 

flu, Zqi] 



(1) 



where R and z are the cylindrical coordinates, p w (0,0) is 
the warm dust density at the center of the Galaxy and h w 
and z w are the exponential scalelength and scaleheight re- 
spectively. The emission of the dust depends on the tem- 
perature T w and on the emission cross-section per unit 
mass of dust n\. In the wavelengths longwards of 60mum, 
the extinction and the emission cross-sections per unit 
mass have the same numerical value because the albedo 
is negligible. Therefore, throughout the paper we will use 
the same symbol k\ to denote the emission and extinc- 
tion cross-sections of the dust. For the temperature of the 
warm dust T w , we assume that it is constant throughout 
the Galaxy and equal to 35 K. This is a reasonable as- 
sumption since the warm dust is mostly associated with 
regions of intense star formation. The spectral energy dis- 
tribution (SED) of a typical star-forming region can be 
approximated by thermal emission from dust with a tem- 
perature of 35 K (see Popescu et al. 120001 and Misiriotis 
et al. I2004|l . For k\ we use the Ry = 3.1 model from 
Weingartner & Draine i200ip . Thus, at a given wavelength 
A, at a given point (R, z) in the Galaxy, the emission per 
unit volume per unit wavelength per steradian -q w (A, R, z) 
can be calculated by 

i] w (X,R,z) = p w (R,z)K X B(X,T w ), (2) 
where B(X,T) is the Planck function. 

3.1.2. The cold dust component 

The density of the cold dust component is described by a 
similar exponential disk according to the following equa- 
tion 



,,,(/?.:) />.(IU))(>xi)(~-M 



(3) 



Here p c (0,0) is the cold dust density at the center of the 
Galaxy and h c and z c are the exponential scalelength and 
scaleheight, respectively. For the temperature of the cold 
dust, we assume that it follows an exponential distribution 
within the Galactic disk, namely, 



T c (R,z) = [T c (0,0)-T oo ]exp 



R_ 



(4) 



Here Iit and zt are the exponential scalelength and scale- 
height for the temperature of the cold dust with T c (0,0) 
being the temperature at the center of the Galaxy and 

Too the temperature of the dust at infinity (in our case 

http://cfa-www.harvard.edu/mmw/MilkyWayinMolClouds.htlmken equal to 3 K). The emission of the cold dust per 
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unit volume per unit wavelength per steradian r] w (X, R, z) 
given by the formula 



is 



7/ c (A, R, z) = p c (R, z)k\B(\, T c (R, z)), 



(5) 



which is identical to Eq. J5J for the warm dust component 
with the dust temperature now varying within the galaxy 
according to Eq. (4). 

Given Eqs.© and (JSJ), the FIR/submm intensity along 
any line of sight /(/, 6, A) can be calculated by 



I(l,b, A) 



Earth 



ds[r) w (\, R, z) + T] c (\, R, z)}, 



(G) 



where I and b are the galactic longitude and latitude, re- 
spectively. The Earth is assumed to lie exactly on the 
galactic plane and at a distance equal to 8 kpc from the 
galactic center (Bahcall & Soneira 1980). The lower limit 
of the integral sq is a point where the line of sight inter- 
sects a "box" inside which we assume that all the galaxy 
is contained. 

In total, the dust spatial distribution 
is described by 10 parameters (namely, 

p w (0,0),h w ,z w ,p c (0,0),hc,z c ,T c (p,0),h T , zt, and 
T w ), which can be constrained using the COBE/DIRBE 
data at 60, 100, 140, and 240 pm. 

3.2. The optical/NIR wavelength range 

Given the dust distribution, as this is inferred by the FIR 
emission, one can calculate the extinction along any line 
of sight at any distance from the Earth. Using this infor- 
mation, we calculate the observed starlight intensity along 
any line of sight at 1.2 and 2.2pm. In particular, the inten- 
sity is given by the following equation (Kylafis & Bahcall 
1987) 



1(1, b, A) 



Earth 



ds (A, i?, z)e 



r 



ds' K\p c {R,z) 



(7) 



-so 



Here n s {\, R, z) is the stellar emissivity at wavelength 
A at position (i?, z) in the galaxy, k\ is the extinction 
cross-section per unit mass at wavelength A (taken from 
Weingartner & Draine 200TJ), and p(R, z) is the dust den- 
sity. To simplify the calculations, we assume that only the 
diffuse cold dust component contributes to the opacity, 
while the warm dust component, being associated mostly 
with the star-forming complexes (see Sect. O , affects the 
opacity locally without affecting the global properties of 
the Galaxy. As we will also see later (Sect.|SJ), the mass of 
the warm dust component is more than 2 orders of magni- 
tude less than that of the cold dust component, making it 
only a very small fraction of the ISM. For the spatial dis- 
tribution of the stellar emissivity we follow Xilouris et al. 
H1997llT§§glll999[l and assume that the stellar distribution 
is described by 



, R \z\ 
■q s {R, z) = rj dlsk exp ( 

ft a Za 



T]bulge 



exp(-7.67 J B 1 / 4 )B 



-7/8 



(8) 



In this expression the first part describes an exponential 
disk and the second part describes the bulge, which in pro- 
jection is the well-known i? 1 / 4 -law fChristensen 11990| . In 
the first term of the expression, r\disk is the stellar emis- 
sivity at the center of the disk, and h s and z s are the 
scalelength and scaleheight, respectively, of the stars in 
the disk. For the bulge, rjbuige is the stellar emissivity per 
unit volume at the center of the bulge, while B is defined 

by 



B 



y/R 2 + z 2 (a/by 

Re 



(9) 



with R e being the effective radius of the bulge and a and 
b being the semi-major and semi-minor axis, respectively, 
of the bulge. 

To model the Galaxy in the NIR wavelengths (and in 
particular at 1.2 and 2.2pm for which we have available full 
maps of the Galaxy), we need to calculate 9 parameters. 
Three of them (p c (0, 0), h c , z c ; the parameters of the cold 
dust component) are constrained from the 60, 100, 140, 
and 240 pm data, as we will see later on (see Sect. 4). The 
remaining 6 parameters (i] dtsk , h s , z s ,i] bu i ge , R e , and a/b) 
are constrained from the 1.2 pm map. Additional modeling 
is also performed on the 2.2 pm map of the Galaxy (see 
Sect.EJ. 

3.3. The molecular and atomic hydrogen 

The two gas phases (HI and H2), being totally indepen- 
dent from the dust and stars, are treated separately from 
the rest of the components. The distribution of the molec- 
ular hydrogen (H2) is assumed to follow the CO emission. 
Although it is rather clumpy, it can be, to a first approx- 
imation, described by an axisymmetric exponential disk 
(e.g., Regan et al. 2001 for the the large-scale distribution 
of CO in other nearby galaxies). Thus, the distribution of 
H2 as a function of the position in the galaxy is given by 



pn 2 (R, z) = pn 2 (0, 0) exp ( - — 

"H 2 Z H 2 



(10) 



where pn 2 (0, 0) is the density of the H2 molecules in the 
center of the galaxy, hn 2 is the scalelength, and zh 2 is the 
scaleheight of the distribution. 

For the atomic hydrogen (HI), we assume a similar dis- 
tribution, but we also introduce an inner truncation radius 
to account for the absence of HI in the central parts of the 
Galaxy, as indicated by the relevant map (see Kalberla et 
al. 12005(1 . The distribution of HI is then given by 



pm{R, z) 



m (0,0)exp -^-ifL , p>R 







, p < Rt 



(11) 



with phi(0, 0) as the HI density at the center of the galaxy 
(of the untruncated disk), R t as the inner truncation ra- 
dius, /ihi as the scalelength, and zbi as the scaleheight of 
the distribution. Although more sophisticated functional 
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Table 1. Masked sources in DIRBE maps. 



Source 


Central b 


Subtracted area 




(degrees) 


(sq. degrees) 


Galactic center 


0, 


10 x 10 


Orion arm 


-80, 


20 x 20 


High latitude 


10, 20 


20 x 20 


Magellanic clouds 


80, -30 


20 x 20 


M31 


-120, -20 


30 x 20 



expressions have been suggested for the description of the 
HI distribution (i.e., van den Bosch 20QOJ); we adopt the 
simplest possible description. 

The model column density of the molecular and the 
atomic hydrogen along any line of sight is then calculated 

by 

pEarth 

I eas (l,b)= / ds Pgas (R,z) (12) 

J s 

with p gas (i?, z) being equal to either the molecular 
Pb_ 2 (R, z) or the atomic pm(R, z) hydrogen distribution 
(see Eqs. (10) and (11), respectively). The molecular hy- 
drogen distribution is then parametrized by three param- 
eters (ph 2 (0j 0), /ih 2 j an( i z h 2 )i while the atomic hydrogen 
distribution, by four (phi(0, 0), /ihi, zjji, and R t ). 

4. The fitting procedure 

Having a computationally intensive radiative transfer 
model with 15 parameters in total to compute (9 for the 
dust distribution and 6 for the stellar distribution; see 
Sect. 3), we proceed in two steps. We first constrain the 9 
parameters that define the dust distribution and the tem- 
perature of the cold dust component by fitting the model 
FIR intensity (Eq.© to the FIR data (i.e., 60, 100, 140, 
and 240/im). Then, having defined the dust distribution, 
we fit Eq. Q to the NIR maps to constrain parameters 
that describe the stellar distribution. 

A third step, treated separately, is to determine the 
parameters that describe the H2 and HI distributions. 

In all three steps, the observed surface brightness (or 
the column density in the case of the gas) is compared with 
the computed surface brightness from the model. Before 
going into a detailed x 2 minimization to find those values 
of the parameters that best describe the Galaxy, we fit 
simple exponential profiles to the surface brightness and 
in directions vertical and parallel to the disk (see Xilouris 
et al. 1997 for more details). In this way, good estimates of 
the geometrical characteristics are derived and are used as 
initial guesses in a % 2 minimization algorithm. The min- 
imization is done using the Steve Moshier C translation 
of the public domain Levenberg-Marquardt solver of the 
Argonne National Laboratories MINPACK mathematical 
library 3 . We always test the uniqueness of the best values 
of the parameters derived by the fit by altering the ini- 
tial values of the fit by as much as ~ 30% — 40%. In all 

3 http://www.netlib.org 



cases the values returned by the fit were the ones that we 
present. 

For all the DIRBE maps, we used the data within the 
latitude range of -40 to 40 degrees. This was done because 
outside this latitude range the signal drops significantly 
with no usable input to the model. Some regions were also 
masked from the DIRBE maps due to their large deviation 
from the Galactic emission (see Table 0. 

5. Results 

Having described the way to perform the fitting of the 
model to the DIRBE maps (Sect.0J, we derive the values 
of the parameters of the model that best fit the data. The 
fitted parameters are presented in Tables [3 El and 0J In 
TableElwe give the values of the parameters that describe 
the dust distribution inside the Galactic disk. In Table 
[3] we present the values of the parameters for the stellar 
distribution as they have been determined from the 1.2 
and 2.2 /jm maps. Finally, in Table^Jwe present the values 
of the parameters for the distribution of the molecular and 
the atomic hydrogen. 

Having determined all the parameters for the dust and 
the stellar distributions, we can now create the model im- 
ages of the Galaxy at each wavelength to compare with 
the observed images. This is what we do in Fig. ^ where 
the observations at different wavelengths (left panels) are 
compared with the model images (right panel). The 1.2, 
2.2, 60, 100, 140, and 240 pm images are shown from top to 
bottom. In these images the x- and y-axes are the Galactic 
longitude and latitude in degrees with the Galactic cen- 
ter at (0,0), while positive and negative longitudes are 
the southern and northern regions of the sky, respectively. 
Positive and negative latitudes are regions above and be- 
low the Galactic plane, respectively. 

As it is obvious from Fig. ^ our model and the ob- 
servations compare quite well for all the different wave- 
lengths. The NIR (1.2, 2.2 //m) maps (two top panels) 
show the characteristic dust lane crossing the bulge along 
the Galactic plane which, at these wavelengths, is poorly 
visible (compared to the optical wavelengths) . The stellar 
disk is well described by the smooth exponential distri- 
bution used in the model (Sect. 0J. The FIR (60, 100, 
140, and 240 pm) maps (third, fourth, fifth, and sixth 
panels from the top, respectively) show the diffuse dust 
emission along the Galactic plane. Given the simplicity of 
the axisymmetric distributions that were used, the model 
quite accurately reproduces the observed emission at these 
wavelengths. 

To get a better view of the goodness of the fit of the 
model to the data, we produce vertical profiles of the 
surface brightness along 10 meridian zones of 10 degrees 
width each. This is shown in Fig. where in each panel 
we overplot the model to the data for each wavelength 
(1.2, 2.2, 60, 100, 140, and 240 /im, as indicated with the 
wavelength value inside each panel) . In each panel (of dif- 
ferent wavelength) the averaged profiles along the longi- 
tude meridians of 0, 5, 10, 20, 30, 60, 90, 120, 140, and 
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Fig. 1. The COBE/DIRBE maps of the Galaxy (left panels) in direct comparison with the fitted model (right panels). 
From top to bottom the 1.2, 2.2, 60, 100, 140, and 240/j,m maps are presented for the Galactic latitude range of -30 to 



30 degrees. The maps were made in Mollweide projection (Snyder 1987). 



180 degrees are presented shifted with each other by a 
factor of 10 in surface brightness for reasons of clarity. As 
one can see, the agreement between the model and the 
observations is quite good. 

We note here that in the first two panels of Fig. |2thc 
absorption by dust is seen as a dip in the first few meridian 
profiles. The dip is more prominent at 1.2 than at 2.2 
yum. The 60 fim emission is mainly due to the warm dust 
that traces regions of star formation. The 240 /im emission 
traces the diffuse cold dust. 

To get a better feeling of the goodness of the fit, we 
present the relative percentage of the residuals as a func- 
tion of the percentage of the Galaxy's area. We do that 
in Table 5 where, for example, for the 1.2 micron band 
we have 29% of the area of the Galaxy with residuals 
less than 10%, 53% of the area with residuals less than 
20%, and 92% of the area of the Galaxy with residuals 
less than 50%. The same holds for the rest of the bands. 
We see that, on average, about 50% (half of the Galaxy's 
area) has residuals less than 20%. Given the large noise 
measurements that exist in the high latitude regions (es- 
pecially in the FIR maps), and also the complexity of the 
real structures, these numbers show how well the model 
fits the real data. 



Using the values of the parameters presented in Table 
13 and adopting the Weingartner & Draine H2()()l(l values 
for the extinction cross-section at 1.2/im and at 2.2/im, 
we can calculate the optical depth (and subsequently the 
extinction) between any two points in the Galaxy. In par- 
ticular, the central face-on optical depth = 2n\z c is 
0.33 at 1.2 /im and 0.17 at 2.2 fim. 

In Fig. |21 we compare the maps of the hydrogen column 
density with the corresponding model. The top panel on 
the left shows the molecular hydrogen map as this is in- 
ferred from the CO observations (Dame et al. 120011 . The 
respective model column density is shown in the top right 
panel. From the comparison between the model and the 
observations it is evident that, despite the large dumpi- 
ness of the data, the large-scale structure is fairly well 
represented by the model. In the bottom left panel we 
show the map of the atomic hydrogen (Kalberla et al. 
2005). The lack of this component in the central part of 
the galaxy shows up as a decrease of the column density. 
The model column density is shown in the bottom right 
panel. In this image, the lack of atomic hydrogen in the 
center of the Galaxy is more prominent due to the sharp 
cut introduced by the model at the truncation radius Rt- 
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Fig. 2. Profiles along Galactic meridian zones (circles) together with the corresponding models (solid lines). Taking 
advantage of the symmetry of the Galaxy, we have folded the maps over the central meridian and extracted 10 meridian 
profiles (centered at 0, 5, 10, 20, 30, 60, 90, 120, 140, and 180 degrees in Galactic longitude). Each profile is averaged 
over a 10 degree zone in Galactic longitude. Due to the folding of the map, the zero longitude profile and the 180 
longitude profile come from the 0-5 and 175-180 meridian zones, respectively. From top left to bottom right we show 
the profiles at 1.2, 2.2, 60, 100, 140, and 240/zm (as indicated in the upper left corner of each panel). 

The comparison between the model and the data is 
shown better in Fig. 0] The left panel shows vertical pro- 
files of the column density of the atomic hydrogen along 10 
meridian zones of width 10 degrees each. The agreement 
between the model and the data is good. The observed 
CO (and subsequently the H2) distribution on the other 
hand is very clumpy as it can be seen from the noise of 
the data in the profiles shown in the right panel of Fig.0] 
The smooth model, however, follows the actual distribu- 
tion fairly well. 

As in the case of the COBE/DIRBE maps, we also 
present the residuals (in percentage) between the model 
and the observations as a function of the Galaxy's area. 
We do that in Table 5 (last two rows). 

Having determined the geometrical characteristics of 
the distributions and the central densities, it is then 
straightforward to compute the total mass of the cold dust 
M c , the warm dust M w , the mass of the molecular hydro- 
gen Mh 2 , and the mass of the atomic hydrogen Mhi by 
integrating the respective three dimensional distributions 



in space. For distributions of the form of Eqs. (1), (3), and 
(10) for the warm dust, the cold dust, and the molecular 
hydrogen, respectively, this integration gives 

M = 47rp(0, 0)zh 2 , (13) 

with h being the scalelength of the distribution, z the 
scaleheight, and p(0, 0) the central density of the material. 
If we also include an inner cut-off in the distribution with 
a truncation radius Rt , as we do for the atomic hydrogen 
(Eq. (11)), the integration then gives 

M = 47rp(0, Q)zh(R t + h) exp (~^\ ■ ( 14 ) 

Using the values for the fitted parameters presented in 
Tables|21and0]and the above equations, we then find M c = 
7.0 x 1O 7 M , M w = 2.2 x 10 5 Af o , A/ H2 = 1.3 x 1O 9 M , 
and A/hi = 8.2 x 10 9 Af Q . The total gas-to-dust mass ratio 
becomes (Mjj 3 + Mm) / (M c + M w ) = 135, which is in fair 
agreement with the value of 160 ± 60 reported in Sodroski 
et al. (dll. 
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Fig. 3. The gas maps of the Galaxy (left panels) in direct comparison with the fitted model (right panels). The upper 
left panel shows the molecular hydrogen distribution (Dame et al. 1200 1[1 . while the atomic hydrogen distribution 
(Kalberla et al. 121) (J5f) is presented in the bottom left panel. Both maps are given for the Galactic latitude range of -30 
to 30 degrees. 
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Fig. 4. Profiles along Galactic meridian zones (circles) together with the corresponding models (solid lines). As in 
the case of the COBE/DIRBE data (Fig. 2), the maps are folded over the central meridian and 10 meridian profiles 
are extracted following the same method described in Fig. [S] The left panel shows the atomic hydrogen profiles, and 
the right panel shows the molecular hydrogen profiles. 

Table 3. Parameters for the stellar distribution (see text for a detailed description of each parameter). 



Parameter 


Unit 


Value 








(1.2/Ltm) 


(2.2^m) 


fjdisk 


erg sec -1 cm~' i Hz _i srad _i 


5.49 x 10" JS 


9.94 x 10 _3S 


lis 


kpc 


2.5 


2.2 


Zs 


kpc 


0.16 


0.12 


fjbulge 


erg sec -1 cm~ 3 Hz -1 srad -1 


2.07 x 10~ 34 


2.03 x 10~ 35 


Re 


kpc 


0.68 


0.79 


a/b 




0.61 


0.63 



6. Discussion 

6.1. Comparison of the model with the FIRAS data 

To further establish the validity of the model, we com- 
pare it with the COBE/FIRAS data. In Fig.EJwe present 
the spectrum of the Galaxy in several directions along 
the Galactic plane averaged over an area of ±10 degrees 
in latitude and 10 degrees in longitude centered on the 
Galactic equator and at 5, 10, 20, 30, 60, 90, 120, and 140 



degrees in longitude. On average, the model is in good 
agreement with the data. The deviation in some cases is 
to be expected due to the non-axisymmetric nature of the 
real distribution of the dust. 

Evidence supporting the previous statement comes 
from a comparison of the "global" SED of the Milky Way 
constructed by averaging the COBE/FIRAS maps in all 
the directions on the sky with the predicted SED from our 
model. This is shown in Fig. El (top panel). 
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Fig. 5. The Spectral Energy Distribution (SED) of the Milky Way averaged between -10 and 10 degrees in Galactic 
latitude and over 10 degrees in Galactic longitude, as indicated on top of each panel. In each panel the filled circles 
are the COBE/FIRAS data and the open circles are the COBE/DIRBE data. The dashed line is the contribution to 
the SED of the warm dust component, the dot-dashed line is the contribution of the cold dust component and the 
solid line is the total modeled SED. 



Table 2. Parameters for the dust distribution (see text 
for a detailed description of each parameter) . 



Parameter 


Units 


Value 


P»(0,0) 


gr cm -3 


1.22 x 10"^' 


h w 


kpc 


3.3 




kpc 


0.09 


Pc(0,0) 


gr cm~ 3 


1.51 x 10" 25 


h c 


kpc 


5.0 


Zc 


kpc 


0.1 


T c (0,0) 


K 


19.2 


flT 


kpc 


48 


ZT 


kpc 


500 



The goodness of the fit can be evaluated by looking 
at the residuals between the model and the observations 
which are presented in the bottom panel of Fig. [S] From 
this plot we see that 82% of the data points on this SED 
show residuals from the modeled SED that are less than 
10%. In particular, 14%, 24%, 47%, 82%, and 93% of the 



Table 4. Parameters for gas distribution (see text for a 
detailed description of each parameter). 



Parameter 


Units 


Value 


Pff 2 (0,0) 


cm" 3 


4.06 


hu 2 


kpc 


2.57 


ZH 2 


kpc 


0.08 


PHl(0,0) 


cm" 3 


0.32 


flHI 


kpc 


18.24 


ZHI 


kpc 


0.52 


Rt 


kpc 


2.75 



data points show residuals from the model which are less 
than 1%, 2%, 5%, 10%, and 20%, respectively. 

This agreement, which does not come from fitting the 
model to the data, but from comparison of the data with 
our model, gives us confidence that our model gives a con- 
sistent representation of the Galaxy in wavelengths ex- 
tending up to 1000 /*m. 

In their analysis, Reach et al. I|1995|) reported emission 
from a very cold component (4-7 K). This component is 
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Table 5. Residuals between the model and the observa- 
tion. For each COBE/DIRBE map we present the per- 
centage of the Galay's area with residuals less than <10% 
, <20%, and <50%. 



Wavelength (jMm) <10% <20% <50% 



1.2 


29% 


53% 


92% 


2.2 


32% 


55% 


87% 


60 


26% 


52% 


88% 


100 


23% 


43% 


81% 


140 


28% 


51% 


81% 


240 


18% 


37% 


75% 


HI 


17% 


38% 


85% 


H 2 


17% 


35% 


67% 



; -180 < Longtitude < 180 




Fig. 6. The global SED of the Milky Way averaged over 
the whole sky (top panel). As in Fig. [S] the filled cir- 
cles are the COBE/FIRAS data and the open circles are 
the COBE/DIRBE data, while the dashed line is the con- 
tribution to the SED of the warm dust component, the 
dot-dashed line is the contribution of the cold dust com- 
ponent, and the solid line is the total modeled SED. In 
the bottom panel we present the residuals (in percentage) 
between model and real measurements. 



not evident in our analysis when looking at the averaged 
SED of the Galaxy (Fig. |5J). However, deviations though 
between the model and the observations do exist when 
looking through various lines-of-sight (Fig.[5j), which can, 
as said earlier, be explained by local deviations of the den- 
sity of the dust material. 



6.2. Temperature distribution 

As already mentioned in Sect. 3.1, the dust in our model 
is described by two components (one with a constant 
"warm" temperature of 35 K, associated with the warm 
HII star-forming regions, and another "cold" temperature 
component). For the distribution of the temperature of 
the cold dust material, we assume a simple exponential 
distribution in both directions, radially and vertically to 
the plane of the Galactic disk (see Eq.0}, with a boundary 
condition that the temperature at large distances is at 3 
K (of course, as noted earlier, the radiative transfer calcu- 
lations only take place within a finite volume where dust 
material exists; see Sect. 3). To better visualize the spatial 
distribution of the temperature of the cold dust compo- 
nent we plot the temperature map seen in Fig.0 This map 
represents a two dimensional slice of the Galactic plane 
that extends up to 20 kpc along and up to 0.3 kpc above 
the Galactic plane with the temperature contours plotted 
every 0.5 K. The radial profiles of the cold dust density 
(in units of 5 x 10 -25 gr/cm 3 ) are also overplotted for 
z=0.0, z=0.1 kpc, and z=0.2 kpc indicating the way that 
the dust material is distributed within this temperature 
field. From this plot we see that the temperature varies 
from 19.2 K at the center to ~ 15 K at the outer parts of 
the Galaxy (~ 15 kpc). This result agrees with the studies 
of Reach et al. (|1995[l . Sodrosky et al. (1994), and Davies 
et al. (1997) showing the similar radial dependence of the 
cold dust component. Concerning the vertical dependence 
of the temperature of the cold component, we see that 
this is more or less constant for different heights above 
the Galactic plane. We note here though that for higher 
opacities this picture could be different, with the dust be- 
ing more effectively heated in places above the disk com- 
pared with the dust in the Galactic plane, where the dust 
is shielded from the diffuse radiation field (see Bianchi et 
al. 2000). 

6.3. Star formation efficiency in the Milky Way disk 
and the Schmidt law 

As already mentioned in the Introduction, a good indica- 
tor of the SFR is the FIR emission (e.g., Hunter et al. 1986; 
Lehnert & Heckman 1996; Meurer et al. 1997; Kennicuut 
1998b; Kewley et al. 2002). To exploit the FIR emission 
we make use of the strong correlation between the SFR 
and the 100 /im flux found in Misiriotis et al. (2004) for 
a large sample of spiral galaxies. Using this relation and 
calculating the total 100 /jm luminosity, we derive a SFR 
of 2.7M yr- 1 for the Milky Way. 

The SFR surface density can then be calculated in the 
disk of the Milky Way, using the same relation presented 
in Misiriotis et al. (2004), and plotted as a function of the 
Galactic radius. This is done in Fig. [SJ where our model 
prediction (solid line) is compared with several estimates 
of the SFR found in the literature (see Boissier & Pratzos 
1999 and references therein). Both the SFR normalized to 
the value at the galactocentric distance of 8 kpc (SFR Q ) 
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Fig. 7. Temperature map of the cold dust component within the Milky Way. Temperature contours (the almost vertical 
black solid lines) are plotted every 0.5 K and highlighted at 15 K and 19 K by a label and a thicker line. The radial 
profiles of the cold dust density are also overplotted for z—0.0 (red solid line), z—0.1 kpc (blue dotted line), and z=0.2 
kpc (green dashed line) in units of 5 x 1CP 25 gr/cm 3 . 



and the star formation rate density are presented in this 
plot as a function of R, with the model showing an excel- 
lent agreement with the existing estimates of the various 
SFR tracers to at least two kpc from the center. The in- 
ner region shows, acknowledging the large uncertainty of 
the few data points available, a decrease of the SFR effi- 
ciency, presumably due to a non-exponential distribution 
of the dust in such small scales. Lack of dust, for exam- 
ple, in this region may be due to the presence of the bar 
structure of the Galaxy. As mentioned earlier though, the 
purpose of this study is to keep the modeling as simple as 
possible with the aim of deriving a valuable description of 
the Galaxy. With this goal in mind we avoid going into a 
more sophisticated modeling of the central part. 

A simple power-law relation between SFR and the gas 
content of external galaxies, introduced by Schmidt ( 1959 ) 
and further explored by Kennicutt (1998a), is well estab- 
lished and tested for large samples of galaxies (Misiriotis 
et al. 12004(1 . This relation is expressed in terms of the SFR 
surface density T,$fr and the gas surface density S gas and 
has the form: 



^SFR - ATt-g. 



(15) 



with A being the absolute SFR efficiency and N the 
slope of the power-law. Kennicutt (1998a) examined two 
samples of galaxies showing normal and starbursting be- 
havior and derived a Schmidt law for normal galaxies 
(with N = 2.47 ± 0.39) and for starburst galaxies (with 
N = 1.40 ±0.13). 

In our study we have, for the first time, the opportunity 
to examine the validity of this law not for a sample of 
external galaxies, by measuring their global properties, 
but for different regions of the same galaxy (the Milky 
Way). This is shown in Fig.[5]with the SFR surface density 
(Esfr) calculated for different radii along the Galactic 
disk, plotted against the gas surface density E gas of the 
same region. These data are shown with the numbers close 



to the points indicating the distance to the Galactic center 
(in kpc). Furthermore, it is evident that for regions within 
the Galaxy with small gas content (at the outer part of the 
disk) the SFR is low, while it gets more intense for regions 
with larger gas density (close to the center of the galaxy) . 
From this plot it is evident that the Galactic Schmidt law 
also follows a power law. A best fit to these data (indicated 
with a solid line in the plot) yields 



J SFR 



Moyr-ipc- 



= 4.45 x 10" 



M Q kpc- 



2.18 



(16) 



The power-law index of 2.18±0.20 that we find is in a very 
good agreement with the value of 2.47 ± 0.39 quoted by 
Kennicutt (1998a) for normal galaxies. This is also evident 
by looking at the distribution of the normal galaxies on 
this plot (Fig. El). The points for the external galaxies 
are well correlated with the measurements for different 
regions within the Milky Way. This indicates the strong 
correlation between the gas content and the SFR activity. 
The striking similarity of the "internal" Schmidt law for 
an individual galaxy (in this case the Milky Way) and the 
"external" Schmidt law for the global properties of normal 
galaxies may indicate a universal gas-SFR correlation not 
only on a large scale (for galaxies as a whole), but also on 
a small scale (for different regions inside a galaxy). 

Another very useful parameter that describes the cur- 
rent star-forming activity of a galaxy is the star formation 
efficiency (e). In Fig.[U]the three parallel, lines correspond 
to constant SFRs per unit gas mass in units of 1%, 10%, 
and 100% per 0.1 Gyr, as described in Kennicutt JT998a). 
It is interesting to notice that the star formation efficiency 
of the Milky Way goes from ^1% per 0.1 Gyr in the outer 
regions (~ 14 kpc) to ~10% per 0.1 Gyr in the center. 
This means that the Galaxy converts ~1% of the gas in 
the outer parts of the disk (~10% in the center) to stars 
over the period of 0.1 Gyr, which roughly corresponds to 
one orbital period of the disk. 
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Fig. 8. The star formation rate (SFR) normalized to the 
value at the galactocentric distance of 8 kpc (SFR©; y-axis 
on the left) and the star formation rate density CEsfr] 
y-axis on the right) as a function of radius in the Milky 
Way. The data correspond to different star formation esti- 
mates and are taken from Lyne et al. (1983J open circles), 
Guesten & Mezger l|1982l triangles), and Guibert et al. 
I|1978l squares; see Fig. 2 of Boissier & Prantzos I1999fl . 
The solid line corresponds to the model prediction. 
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Fig. 9. Star formation rate density as a function of gas 
surface density in the Milky Way disk. The solid circles 
indicate the measurements for the Milky Way (for dis- 
tances along the disk as indicated by the numbers close 
to the solid circles), while the solid line is the best fit of 
a Schmidt law to these points (see the text for details). 
The open circles are the global measurements for a sam- 
ple of external galaxies presented in Kennicutt {1998a). 
The three parallel dashed and dotted lines correspond to 
constant star formation efficiencies of 1%, 10%, and 100% 
per 0.1 Gyr. 



The average global efficiency of star formation in the 
Milky Way can be calculated by dividing the SFR over 
the period of 0.1 Gyr (2.7 x 10 8 M Q ; see beginning of Sect. 
6.3) by the total gas mass (M H2 + M HI = 9.5 x 10 9 ). 
This results in 2.8% per 0.1 Gyr or, in other words, that 
our Galaxy spends 2.8% of its gas to create stars over 0.1 
Gyr. This is at the lower end of the median star formation 
efficiency of typical present-day spiral galaxies, which is 
4.8% (Kennicutt QMS). 

6.4. Comparing the model with optical data 

To check the validity of our model in the optical 
wavelengths, we use the Database of Galactic Classical 
Cepheids (Fernie et el. I1995J) . This database contains in- 
formation about over 500 classical Cepheids in the Galaxy 
with determined distances ranging from a couple of hun- 
dred parsec up to ~ 6.5 kpc. The E(B-V) color excess 
reported in this catalog for each star was then checked 
against the E(B-V) determined with our model by cal- 
culating for each Cepheid the extinction in the B- and 
V-bands. The B and V extinction maps for the Galaxy 
were created by assuming the values of h c and z c given 
in Table 2 and using the values of the extinction coeffi- 
cient for the B- and V-bands calculated in Weingartner & 
Draine ( 2001). The comparison between the observed data 
and the model predictions for E(B-V) is shown in Fig. 1101 
with the x-axis being the B-V color excess observed and 
the y-axis the B-V color excess predicted by the model. 
From this plot it is evident that the correlation between 
the observed and the predicted measurements is extremely 
good with the average of the points, binned every 0.2 mag- 
nitudes of E(B-V) on both axes, falling, within the error 
bars, along the diagonal of the plot. This is also evident by 
looking at the individual observations that are scattered 
symmetrically with respect to the diagonal. 

The scatter of the points indicates that in many cases 
our model either overestimates or underestimates the 
amount of dust along the line of sight. This is clear ev- 
idence of the clumpy distribution of the dust. However, 
the fact that the opacity overestimated lines of sight are 
about as many as the underestimated ones ensures that on 
average our smooth model is consistent with the observed 
extinctions. 

7. Conclusions 

In the present study a consistent model of the ISM for a 
wide range of the spectrum (from optical wavelengths to 
millimeter wavelengths) is presented. Using a simple ax- 
isymmetric three-dimensional radiative transfer model for 
the dust and the stars in the Galactic disk, we derive the 
set of the parameters that best describes the FIR emission 
observed by COBE/DIRBE at 60, 100, 140, and 240 ^m. 
Having constrained the dust distribution, we then fit the 
NIR COBE/DIRBE observations at 1.2 and 2.2 ^m and 
determine the NIR properties of the Galaxy. The scale- 
length for the stars at 1.2 and 2.2 ^m is found to be 2.5 and 
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Fig. 10. The B-V color excess for the sample of the 
Galactic Classical Cepheids (Fernie et el. 1995; x-axis) 
compared with the color excess along the same directions 
as derived from the model. 



2.2 kpc respectively. This comes in agreement with find- 
ings of recent studies (see Drimmel & Spergel 2001 and ref- 
erences therein). Having calculated both the NIR bands, 
we see the trend of increasing scalelength with decreas- 
ing wavelength. This trend is also found in other sudies of 
nearby galaxies (e.g., De Jong 1996; Xilouris et al. 1999). 
We find that the ratio of the stellar to the dust scaleheight 
is ~ 1.7, while the scalelength of the dust is about twice 
that of the stars (as inferred from the NIR observations). 
Given the fact that we expect larger scalelength values 
in the optical (by as much as 20%), this scalelength ra- 
tio could become ~ 1.4 — 1.5 in the optical wavelengths 
(see Davies et al. 1997; Xilouris et al. 1999). In addition 
to the COBE data, we also model the 3D distribution of 
the atomic and molecular hydrogen as they reveal them- 
selves from the 21 cm and the CO emissions, respectively. 
We then compare our model predictions (for wavelengths 
up to 1000 jum) with the observations by COBE/FIRAS. 
The global SED of the Galaxy (averaged over the entire 
Galaxy) is in excellent agreement with the model, while 
the model shows small deviations from the SED of specific 
directions on the sky. These could be explained by assum- 
ing local deviations of the density of the dust (dumpi- 
ness of the dust). The star formation rate as a function 
of Galactocentric distance is then derived by the model 
and is found to be in excellent agreement with various 
star formation tracers within the Galaxy. Comparing the 
SFR surface density with the gas surface density, an "in- 
ternal" Schmidt law is derived for the regions along the 
Galactic disk that is almost identical to the "external" 



Schmidt law derived by Kennicutt ( 1998a) for typical spi- 
ral galaxies. This could be an indication of a universal gas- 
SFR behavior not only on large scales (for the galaxies as 
a whole), but also on small scales (within an individual 
galaxy). The star formation efficiency of the Milky Way is 
found to decrease from 10% to 1% per 0.1 Gyr when go- 
ing from the center to the outer parts of the Galactic disk. 
Finally the model predictions in the optical wavelengths 
(B- and V-bands) are compared with existing observations 
of Cepheid stars (in terms of their observed E(B-V) color 
excess) and a good correlation (within the statistical er- 
rors) is found. 
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